Interfacial confinement in core-shell nanowires due to high dielectric mismatch 
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We theoretically investigate the role of the dielectric mismatch between materials on the energy 
levels and recombination energies of a core-shell nanowire. Our results demonstrate that when 
the dielectric constant of the core material is lower than that of the shell material, the self-image 
potential pushes the charge carriers towards the core-shell interface, in such a way that the ideal 
confinement model is no longer suitable. The effects of this interfacial confinement on the electronic 
properties of such wires, as well as on its response to applied magnetic fields, are discussed. 

PACS numbers: 



Great attention has been devoted to the investiga- 
tion of the electronics and optical properties of core-shell 
nanowires (NW). In particular, the applications of these 
low dimensionality structure in optoelectronic and pho- 
tonic devices are of interest for the electronics industry, 
and much effort has been dedicated to their fabrication. 

In addition, studies of low dimensional systems sur- 
rounded by high dielectric constant materials continue to 
attract attention from many researchers 0, 0] towards a 
continuation of the Moore's law. Recently, wire diame- 
ters of a few nanometers were experimentally achieved, 
[3] and carrier confinement effects in these nanowires have 
been reported with different levels of sophistication. [§- 

This work aims to investigate the dielectric mismatch 
effects on core-shell NW, focusing on the possibility 
of interfacial confinement of the carriers. As for the 
model structure, we consider a semiconductor cylindri- 
cal nanowire (core region) of radius R, surrounded by a 
different material (shell region). The interface between 
core and shell regions is assumed to be abrupt, i.e. the 
materials parameters change abruptly from the core to 
the shell regions. For the heterostructure materials con- 
sidered in this paper, the bands mismatch creates a high 
potential barrier for the charge carriers in the shell, lead- 
ing to a short penetration of the wave functions in this 
region. This rules out the role of the shell width on 
the energy states of the NW, since the wave functions 
does not reach the outer edge of the shell. The nanowire 
electronic structure is obtained by solving numerically a 
Schrodinger-like equation within the adiabatic approach 
and the effective mass framework. [ll[ The total con- 
finement potential Vj,(fti) — AEi(pi) + Sj(/9j) is given 
by the sum of band edges discontinuities A£^(/dj) and 
the self-energy potential Sj(pi), where i = e,lh,hh rep- 
resents the carrier types (electron, light hole and heavy 
hole, respectively). The latter term appears due to the 
dielectric mismatch, and its calculation is based on the 
method of the image charges. The details of this calcu- 
lation can be seen in our supplementary material. 11 1 



is given by Eq.(19) of Ref. |l0|, whereas when the carrier 
is located in the shell region, the self-energy potential 
inside this region can be obtained using a similar expres- 
sion, just by changing the modified Bessel functions of 
the first type, I^{kp e ) and I^(kph), by those of the sec- 
ond type, K^kpe) and K^(kph), respectively. 

The electron-hole recombination energy, E R ~ h =E G + 
E e + Eh (h = Ih or hh), given by the sum of the 
band gap energy Eq on the core region with the elec- 
tron E e and hole Eh confinement energies, is calculated 
for the radial ground state quantum number n = 1 
and zero angular momentum I — 0, allowing us to 
analyze the effects caused by the self-energy potential 
on the ground state energy of the electron-hole pair. 
Figure QJa) shows the role of the dielectric mismatch 
on the value of the electron-hole recombination energy 
Eji^.rpA-Q, as a function of the 
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ratio e r = ei/s-2 between dielectric constants of the core 
(e%) and shell (£2) materials, considering several values 
of nanowire radii R, for different materials. This quan- 
tity denotes the difference between the recombination en- 
ergies with {E R ^.( Pi) ^) and without (-E fljE . (p4 ) =0 ) the 
self-energy corrections. The analysis of this difference 
gives us the advantage of excluding the effect of quan- 
tum confinement due to the band edge mismatch AEi, 
remaining exclusively the confinement due to the dielec- 
tric mismatch. The material parameters are the same as 
in Ref. [H 

As in the case of quantum wells, shown by Pereira et 



12|, if the dielectric constant in core region e± is 



In a nutshell, the self-energy potential T,i(pi) inside the 
core region, due to a carrier located in the core region, 



smaller than that of the shell region £2, the potential 
in the core region is attractive and its contribution to 
the recombination energy is negative AEr < 0. On the 
other hand, for a larger dielectric constant in core region, 
as compared to the one in the shell region, the potential 
in the core region is repulsive, so that its contribution 
to the recombination energy is positive AEr > 0, as 
clearly shown in Fig. [TJa). The resulting AE% for 
both light and heavy holes are the same, as this quantity 
only contains the effects of Tii(pi) on the carrier charge. 
On the other hand, by plotting E e R ~ h , it is possible to 
observe the difference between e — Ih and e — hh pairs, 
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FIG. 1: (a) Effect of the self-energy potential on the recombi- 
nation energy of an electron-hole pair, for different wire radii: 
4 nm (x), 6 nm (•), 8 nm (a) and 10 nm (■). (b) Ground 
state energies as a function of the core radius for carriers con- 
fined in a GaN/Hf02 core-shell heterostructure. 
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FIG. 2: (a) Overlap between the electron (el) and the light 
(Ih) and heavy (hh) hole wave functions as a function of 
the core radius R, for several heterostructure materials with 
different ratios between dielectric constants e r , represented 
by different colors, (b) Effective confinement potentials 
(solid) for R = 30 nm in the GaN/HfO^ case for electrons 
and heavy holes, along with some examples of their wave 
functions. Electron wave functions are shown with (dashed) 
and without (dashed-dotted) taking image charge effects into 
account, for comparison. Hole wave functions (for £ 7^ 0) 
are depicted for R = 15 nm (blue dashed), 20 nm (red dash 
dotted), 25 nm (green dotted) and 30 nm (black short dash 
dotted). 



caused by the different light hole and heavy hole effective 
masses. 

Notice that the AlGaAs/GaAs heterostructure has 
e r ~ 1, i.e. the dielectric mismatch in this case is neg- 
ligible. In fact, one can verify in Fig. (Ua) that for 
AlGaAs/GaAs, AEr » for any value of wire radius. 



AlGaAs/GaAs core-shell wires have been widely exper- 
imentally studied for many reasons, e.g. the possibility 
of obtaining lattice matched defect-free samples, where 
the AlGaAs shell reduces the surface-related nonradia- 



tive recombination, by passivating the GaAs core. [13 1 
However, the results in Fig. 1(a) suggest that this het- 
erostructure is actually not the best choice for investigat- 
ing dielectric mismatch effects. As we will demonstrate 
onwards, the case of GaN/Hf02 heterostructures is par- 
ticularly interesting. 

The growth of nitride-based cylindrical core wires, Q 
as well as Hf02 shells, @ has been reported recently, 
which makes us believe that such a structure could be 
experimentally realized in near future. Figure[IJb) shows 
the ground state energies E (i.e., for n = 1 and I = 0) 
for each carrier in a core-shell GaN/HfO^ heterostruc- 
ture as a function of the core radius. As usual, as the ra- 
dius increases, the carriers confinement becomes weaker, 
leading to a reduction in the ground state energy, which 
is clearly observed for the electron and the light hole 
in this case. Conversely, the heavy hole energy shows 
a non-monotonic dependence on the radius, namely, it 
decreases for smaller radii but starts to increase as the 
radius becomes larger, suggesting a different confinement 
regime for this carrier. In fact, such a behavior is analo- 
gous to the one reported in the literature for GaN/Hf02 
quantum wells. 1J| In the quantum wells case, such a 
non-monotonic behavior of the confinement energy was 
demonstrated to be a consequence of the fact that e r < 1 
in this heterostructure, leading to an attractive poten- 
tial which pulls the carrier towards the interface between 
the materials. Analogously, this suggests the existence of 
an interfacial confinement of the heavy hole in the NW. 
For smaller radii, the quantum confinement due to the 
bands mismatch is still dominant, but as the radius is 
enlarged, the energy decreases and eventually enters the 
negative energy domain, where the hh becomes bound 
to the core-shell interface due to the self-image potential. 
Being compressed towards the interface, the hh experi- 
ences an increase in its energy as the radius is further 
enlarged. 

Notice that the sign of the confinement energy indi- 
cates whether a charge carrier is confined in the core or 
at the interfacial region. Without dielectric mismatch, 
or even for e r > 1 , the confinement energy is always pos- 
itive, since the bottom of the confinement potential Vj. 
is always either zero or larger. Conversely, the presence 
of a e r < 1 mismatch is responsible for negative cusps on 
the potential at the interfacial region. EH Hence, a 
negative confinement energy suggests an interfacial con- 
finement. In fact, the heavy hole is the only carrier with 
negative energy in Fig. 1 (b) and it is indeed the only 
one exhibiting non-monotonic behavior of the confine- 
ment energy as the radius increases. 

The consequences of the interfacial confinement are nu- 
merous. For instance, for some combinations of core-shell 
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materials, one might find holes in the interface and elec- 
trons in the core. Such an electron-hole spatial separa- 
tion leads to a reduction of the overlap between their 
wavefunctions, reducing the oscillator strength and, con- 
sequently, the recombination rate. The overlaps between 
electron and light (Ih, top) and heavy (hh, bottom) hole 
wave functions are shown as a function of the core wire 
radius R in Fig. 2 (a), for the same combinations of 
core and shell materials in Fig. 1 (a). For all the het- 
erostructures shown, the overlap value is slightly lower 
than unit for small radii. This is just a consequence of 
the fact that wave functions for different charge carri- 
ers penetrate with different depths into the barriers, and 
this effect is more pronounced when the NW radius is 
small. In the AlGaAs/GaAs (blue), GaAs/ZnSe (green), 
and Si/Si02 (black) cases, where the dielectric constants 
ratios are e r > 1, the overlaps for both el — Ih and el — hh 
pairs simply converge to 1 as the radius increases, which 
means that electrons and holes are practically equally 
distributed in space for larger wire radius. However, the 
overlaps in the e r < 1 heterostructures exhibit a max- 
imum at moderate R (around 10-20 ran), and decrease 
after this value, indicating different distributions of elec- 
tron and holes functions for large R. Although this ef- 
fect is also present for the Si/Hf02 heterostructure, it is 
much stronger in the GaN/Hf02 case (red), specially for 
the el — hh pair. 

In order to help us to understand the heavy hole inter- 
facial confinement and the consequently different spatial 
distribution of charge carriers, Fig. [2] (b) illustrates the 
the total confinement potential (black solid), along with 
the normalized ground state wave functions, for elec- 
trons (top) and heavy holes (bottom) in a R = 30 nm 
GaN/Hf02 core-shell NW. The electron wave function 
considering dielectric mismatch effects (black dashed) is 
still confined mostly in the vicinity of the central axis of 
the core, even at such a large radius, and its width is 
wider than that without dielectric mismatch (blue dash 
dotted). On the other hand, the heavy hole wave func- 
tions, depicted in the bottom panel of Fig. [2](b) for differ- 
ent values of the core wire radius R, are clearly confined 
at the interface, due to the cusp formed in this region 
by the self-energy potential (see black solid line). No- 
tice that as the radius increases, the hh wave function is 
squeezed towards the interface, explaining the increase of 
the energy in the heavy hole Eq curve in Fig. 1 (b). Due 
to its larger effective mass, as compared to that of the 
other charge carriers, the heavy hole is indeed expected to 
be more strongly confined in the interfacial cusps. More- 
over, heavy holes in GaN/Hf02 also have larger effective 
mass as compared to any charge carrier in Si/Hf02, ex- 
plaining why the interfacial confinement effect is weaker 
in the latter case. 

Another interesting effect arising from the interfacial 
confinement in quantum wires comes from their topology: 
systems where the charge carriers are confined around a 
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FIG. 3: Confinement energies as a function of the magnetic 
field for (a) electrons, (b) light and (c) heavy holes, consider- 
ing n — 1 and different values of angular momentum index I, 
in a R = 30 nm GaN/Hf0 2 NW. (d) The results for a heavy 
hole in the same system, but neglecting image charge effects. 



core are known to produce an interesting effect when a 
magnetic field is applied perpendicular to the confine- 
ment plane, namely, angular momentum transitions oc- 
cur, even for the ground state, as the magnetic field in- 
tensity increases, which is reminiscent of the Aharonov- 
Bohm (AB) effect. |9| In order to illustrate this, the 
electron and hole energy behaviors under such a mag- 
netic field were numerically [9j calculated and shown in 
Fig. |] as a function of the field amplitude B for the 
GaN/Hf0 2 NW. All the carriers have E < 0, suggest- 
ing that their ground states are all interfacially confined. 
A large (R — 30 nm) radius is considered, in order to 
enhance the interfacial confinement effect (see Fig. 2 
(b)) and reduce the AB period, since this period is in- 
versely proportional to the average radius of the carrier 
wave function. Indeed, the ground state energy exhibits 
AB oscillations for all carriers in this case, although they 
are much more evident for the hh. Due to their consid- 
erably smaller effective masses, fl5| el and Ih are just 
weakly confined in the interface, and still have a large 
wave function tail spreading inside the core, explaining 
the large AB period. The first ground state transition 
for these carriers is not so visible in Figs. [3] (a, b), and 
occurs at B w 6 T. On the other hand, for the hh, the 
AB transitions are quite clear, occurring at smaller mag- 
netic fields, with a period B w 2.5 T. One easily verifies 
that the results neglecting the self-energy term in Fig. 
|3] (d) are qualitatively different, showing a more conven- 
tional behavior with the magnetic field for a core-shell 
NW, namely, without angular momentum transitions for 
the ground state. 

The interfacial confinement is also expected to have an 
important role on the electrons mobility along the core 
axis, since this property depends on the electrons density 
in the central region of the wire, which is suppressed 
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in the case of interfacial confinement. However, a more 
detailed investigation on this issue is left for future works. 

In summary, we have theoretically investigated the 
charge carriers confinement in core-shell nanowires with 
strong dielectric mismatch. Our results predict that, for 
specific configurations of the system, the carriers may 
be confined at the core-shell interface. Such interfacial 
confinement leads to drastic modifications on the elec- 
tronic properties of the NW, specially under an applied 
magnetic field, where angular momentum transitions oc- 
cur for the ground state, due to the AB effect. A de- 
crease in the oscillator strength of the electron-hole pairs 
in e r < 1 core-shell quantum wires is also predicted for 
larger wire radii, which directly affects their recombina- 
tion rates. We believe that our results will spur on future 
experimental investigations on core-shell wires made out 
of high-k materials, contributing for a better understand- 
ing of these systems. 
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SUPPLEMENTARY MATERIAL 



Our system consists of a semiconductor cylindrical 
nanowire (core region) of radius R, surrounded by a dif- 
ferent material (shell region). The materials parameters 
are assumed to change abruptly from the core to the shell 
regions. Within the adiabatic approach and the effective 
mass framework, the Hamiltonian for each carrier i = el 
(electron), hh (heavy hole) or Ih (light hole) is given by 





Hi = ( Pi 



Pi - qA) + Vj (>i) ■ 



(1) 



where Pi (= —iKV) is the momentum operator and rtii is 
the effective mass for the i-th carrier, which depends on 
the radial coordinate, due to the different masses in the 
core and shell materials. The wire is assumed to be in- 
finitely long in the z-direction, allowing us to use a plane 
wave solution exp(ik z Zi) for this direction, leading to an 
energy E z — h 2 k1/2m,i for this direction, which is set as 
the energy referential in our calculations, i.e. E z = 0. 
Besides, we take advantage of the circular symmetry of 
our system and choose exp(U9i) as the solution for the 
angular coordinate, where I is the angular momentum in- 
dex. In the presence of a B = Bz magnetic field, taking 
the symmetric gauge A = Bpi9i/2 for the vector poten- 
tial, the Hamiltonian eventually assumes the form 



Hi = - — 



1 d 
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with 



C m {kR, e r ) 



K m (kR) K' m (kR) 



I m (kR) K' m (kR) - e r I' m (kR) K m (kR) : 



(5) 



where I m (x) (K m (x)) is the m-th order modified Bessel 
function of the first (second) kind and e r = £\/e2- 

When the carrier is located in the shell, the Poisson 
equation to be solved in this region is 



e 2 W(r,r') = -e5(r-r'), 



(6) 



where £2 is the dielectric constant of the shell region, and 
the solution V ou t (r, r') for this equation is given by 



Vout (r,r') = 



\r — r 7r 



m— — 00 



x / dkcos[k(z ~ z')]D m (kR,e r ) K m (kp)K m (kp') 




(7) 



with 



The total confinement potential Vj*(pi) — AEi(pi) + 
Tii(pi) is given by the sum of band edges discontinuities 
AEi(pi) and the self-energy potential Sj(pj). The lat- 
ter term appears due to the dielectric mismatch, and its 
calculation is based on the method of the image charges, 
similarly to Refs. Q, 0. The self-energy potential Si(pi) 
inside the core region, due to a carrier located in the core 
region, is calculated by solving the Poisson equation 



D m (kR,e r ) 



I m (kR)I' m (kR) 



I m (kR)K' m (kR) - e r I' m (kR)K m (kR) 



(8) 



The self-energy potential comes from all the interac- 
tions between the carrier and its own image charges, 
which in the core region is calculated by solving the in- 
tegral 



£iV 2 t/(r,r') = -eS (r - r') , 



(3) 



V?(Pi) = g / Hr-n)Vi n (r,ri)dr, (9) 



where £\ is the dielectric constant of the core region, and 
the solution V in (r, r') in this region is given by 



V in (r,r') = 



Awe 1 



„ m=+oo 

1 ^ + -(£,-1) e ™ 

r — r 7r z — ' 

1 1 m— — 00 



(e-e>) 



that results on 



e z t2? r 

Sj" (ft) = ( £ r - 1) E / dkC '™ ( fcjR > £ '') ^Pi) 



(10) 



2 



for electrons and holes (i = el, Ih, hh) in this region. In 
the shell region, the self-energy potential comes from the 
integral 

s° Mt (Pi) = | / S(r-n) Vout (r, n) dr, (11) 
that results on 



by Sf(pi) in Eq. (HUJ) and Ef^OO in Eq. @3) are 
eventually obtained by numerically solving the integrals, 
and add on V^(pi). Finally, the Schrodinger equation 
for the Hamiltonian Eq. ([2]) is discretized in a finite dif- 
ferences scheme, Q and the resulting matrix eigenvalues 
problem is numerically solved by standard computational 
routines, in order to obtain the energy spectrum of the 
quantum wire system. 



^i Ut {pi) = fa- - 1) E / dkD ™ (kR,e r )Kl(k Pl ) 

m— — oo q 

(12) 

Notice the solution in the first case can be compared 
to the second case just by changing the modified Bessel 
functions of the first type, I^ikpi) by those of the sec- 
ond type, K^kpi). The self-energy potentials defined 
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